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Abstract 



In this chapter, we propose a new practical codification of the elements 
(/3 , of the Venn diagram in order to easily manipulate the focal elements. In 

O ■ order to reduce the complexity, the eventual constraints must be inte- 

grated in the codification at the beginning. Hence, we only consider a 
reduced hyper power set Df that can be 2 or D . We describe all 
^ ' the steps of a general belief function framework. The step of decision is 

fT^ ' particularly studied, indeed, when we can decide on intersections of the 

00 . singletons of the discernment space no actual decision functions are easily 

^T ' to use. Hence, two approaches are proposed, an extension of previous 

C^ I one and an approach based on the specificity of the elements on which to 

decide. 
f^ ' The principal goal of this chapter is to provide practical codes of a 

OO , general belief function framework for the researchers and users needing 

^^ ' the belief function theory. 

Keywords: DSmT, practical codification, DSmT decision, 
low complexity. 

X 

. P. , 1 Introduction 

Today the belief function theory initiated by [6, 26] is recognized to propose 
one of the more complete theory for human reasoning under uncertainty, and 
have been appUed in many kinds of apphcations [32] . This theory is based on 
the use of functions defined on the power set 2® (the set of all the subsets 
of O), where Q is the set of considered elements (called discernment space), 
whereas the probabilities are defined only on 0. A mass function or basic belief 
assignment, m is defined by the mapping of the power set 2 onto [0, 1] with: 



E '"W = 1- (1) 



xe2'- 
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One element X of 2®, such as ■m{X) > 0, is called focal element. The set of 
focal elements for m is noted J-m- A mass function where O is a focal element, 
is called a non-dogmatic mass functions. 

One of the main goal of this theory is the combination of information given 
by many experts. When this information can be written as a mass function, 
many combination rules can be used [23]. The first combination rule proposed 
by Dempster and Shafer is the normalized conjunctive combination rule given 
for two basic belief assignments mi and m2 and for all X G 2®, X 7^ by: 

"^Ds(^) = Y3T E ^M)^2{B), (2) 

AnB=x 

where k — > mi{A)in2{B) is the inconsistence of the combination. 

AnB=0 

However the high computational complexity, especially compared to the 
probability theory, remains a problem for more industrial uses. Of course, higher 
the cardinality of is, higher the complexity becomes [38]. The combination 
rule of Dempster and Shafer is #P-complete [25]. Moreover, when combining 
with this combination rule, non-dogmatic mass functions, the number of focal 
elements can not decrease. 

Hence, we can distinguish two kinds of approaches to reduce the complexity 
of the belief function framework. First we can try to find optimal algorithms in 
order to code the belief functions and the combination rules based on Mobius 
transform [18, 33] or based on local computations [28] or to adapt the algorithms 
to particulars mass functions [27, 3]. Second we can try to reduce the number 
of focal elements by approximating the mass functions [37, 36, 4, 9, 16, 17], that 
could be particularly important for dynamic fusion. 

In practical applications the mass functions contain at first only few focal 
elements [7, 1]. Hence it seems interesting to only work with the focal elements 
and not with the entire space 2^^ . That is not the case in all general developed 
algorithms [18, 33]. 

Now if we consider the extension of the belief function theory proposed by 
[10], the mass function are defined on the extension of the power set into the 
hyper power set D^ (that is the set of all the disjunctions and conjunctions 
of the elements of O). This extension can be seen as a generalization of the 
classical approach (and it is also called DSmT for Dczcrt and Smarandache 
Theory [29, 30]). This extension is justified in some applications such as in 
[20, 21]. Try to generate D^ is not easy and becomes untractable for more than 
6 elements in O [11]. 

In [12], a first proposition have been proposed to order elements of hyper 
power set for matrix calculus such as [18, 33] made in 2®. But as we said herein, 
in real applications it is better to only manipulate the focal elements. Hence, 
some authors propose algorithms considering only the focal elements [9, 15, 22]. 
In the previous volume [30], [15] have proposed Matlab^ codes for DSmT hybrid 
rule. These codes are a preliminary work, but first it is really not optimized for 
Matlab and second have been developed for a dynamic fusion. 

Matlab is certainly not the best program language to reduce the speed of 
processing, however most of people using belief functions do it with Matlab. 



Matlab is a trademark of The MathWorks, Inc. 



In this chapter, we propose a codification of the focal elements based on a 
codification of O in order to program easily in Matlab a general belief function 
framework working for belief functions defined on 2® but also on D^ . 

Hence, in the following section we recall a short background of belief function 
theory. In section 3 we introduce our practical codification for a general belief 
function framework. In this section, we describe all the steps to fuse basic belief 
assignments in the order of necessity: the codification of Q, the addition of the 
constraints, the codification of focal elements, the step of combination, the step 
of decision, if necessary the generation of a new power set: the reduced hyper 
power set Df and for the display, the decoding. We particularly investigate the 
step of the decision for the DSmT. In section 5 we give the major part of the 
Matlab codes of this framework. 



2 Short background of belief functions theory 

In the DSmT, the mass functions m are defined by the mapping of the hyper 
power set D^ onto [0, 1] with: 

E ^{X)^l, (3) 

with less terms in the sum than in the equation (3). 

In the more general model, we can add constraints on some elements of 
D®, that means that some elements can never be focal elements. Hence, if we 
add the constraints that all the intersections of elements of Q arc impossible 
(i.e. empty) we recover 2®. So, the constraints given by the application can 
drastically reduce the number of possible focal elements and so the complexity 
of the framework. On the contrary of the suggestion given by the fiowchart on 
the cover of the book [29] and the proposed codes in [15], we think that the 
constraints must be integrated directly in the codification of the focal elements 
of the mass functions as we shown in section 3. Hereunder, the hyper power 
set D taking into account the constraints is called the reduced hyper power set 
and noted Df. Hence, Df can be I?®, 2®, have a cardinality between these 
two power sets or inferior to these two power sets. So the normality condition 
is given by: 

E HX) = l. (4) 

Once defined the mass functions coming from numerous sources, many com- 
bination rules are possible (see [5, 31, 20, 35, 23] for recent reviews of the com- 
bination rules). The most of the combination rules are based on the conjunctive 
combination rule, given for mass functions defined on 2® by: 

s 

"^c(X)= E W^AYi)^ (5) 

Yin...r\Y^=x ]=i 

where Yj £ 2® is the response of the source j, and mj{Yj) the corresponding 
basic belief assignment. This rule is commutative, associative, not idempotent, 
and the major problem that try to resolve the majority of the rules is the 
increasing of the belief on the empty set with the number of sources and the 



cardinality of Q [19]. Now, in D® without any constraint, there is no empty 
set, and the conjunctive rule given by the equation (5) for all X G D^ with 
Yj S Df can be used. If we have some constraints, we must to transfer the 
belief mc(0) on other elements of the reduced hyper power set. There is no 
optimal combination rule, and wc cannot achieve this optimality for general 
applications. 

The last step in a general framework for information fusion system is the 
decision step. The decision is also a difficult task because no measures are 
able to provide the best decision in all the cases. Generally, we consider the 
maximum of one of the three functions: credibility, plausibility, and pignistic 
probability. Note that other decision functions have been proposed [13]. 

In the context of the DSmT the corresponding generalized functions have 
been proposed [14, 29]. The generalized credibility Bel is defined by: 

Be\{X)= Y. ^(^) (6) 

YeDf,YCX,Y^(D 

The generalized plausibility PI is defined by: 

pi(x)= Y. ^(^) (^) 

YeDf,xnY^(/} 

The generalized pignistic probability is given for all X E Df, with X 7^ is 
defined by: 

GPT(X)= J2 ^^f^MY), (8) 

where Cm (X) is the DSm cardinality corresponding to the number of parts of X 
in the Venn diagram of the problem [14, 29]. Generally in 2®, the maximum of 
these functions is taken on the elements in 0. In this case, with the goal to re- 
duce the complexity we only have to calculate these functions on the singletons. 
However, first, there exist methods providing decision on 2® such as in [2] and 
that can be interesting in some application [24], and secondly, the singletons 
are not the more precise elements on Df . Hence, to calculate these functions 
on the entire reduced hyper power set could be necessary, but the complexity 
could not be inferior to the complexity of D^ and that can be a real problem if 
there are few constraints. 



3 A general belief function framework 

We introduce here a practical codification in order to consider all the previous 
remarks to reduce the complexity: 

• only manipulate focal elements, 

• add constraints on the focal elements before combination, and so work on 

• a codification easy for union and intersection operations with programs 
such as Matlab. 



We first give the simple idea of the practical codification for enumerating the 
distinct parts of the Venn diagram and so a codification of the discernment space 
0. Then we explain how simply add the constraints on the distinct elements 
of O and so the codification of the focal elements. The subsections 3.4 and 
3.5 show how to combine and decide with this practical codification, giving a 
particular refiexion on the decision in DSmT. The subsection 3.6 presents the 
generation of Df and the subsection 3.7 the decoding. 

3.1 A practical codification 

The simple idea of the practical codification is based on the affectation of an 
integer number in [1;2" — 1] to each distinct part of the Venn diagram that 
contains 2" — 1 distinct parts with n = |8|. The figures 1 and 2 illustrate 
the codification for respectively 6 — {6*1,^27^3} and 9 = {01,92,0^,64} with 
the code given in section 5. Of course other repartitions of these integers are 
possible. 




Figure 1: Codification for 6 = {Oi,02,0'i}. 

Hence, for example the element Oi is given by the concatenation of 1, 2, 3 
and 5 for \Q\ — 3 and by the concatenation of 1, 2, 3, 4, 6, 7, 9 and 12 for 
|e| = 4. We will note respectively 6*1 = [1 2 3 5] and 6*1 = [1 2 3 4 6 7 9 12] for 
|8| = 3 and for |8| — 4, with increasing order of the integers. Hence, 8 is given 
respectively for |8| =3 and |8| = 4 by: 

8 = {[1 2 3 5],[1 2 4 6],[1 3 4 7]} 

and 

e = {[1 2 3 4 6 7 9 12], [1 2 3 5 6 8 10 13], [1 2 4 5 7 8 11 14], [1 3 4 5 9 10 11 15]}. 

The number of integers for the codification of one element ^i G 6 is given by: 

n-l 

l + Y.^n-1, (9) 




Figure 2: Codification for 6 = {61,62, 63,64}. 



with n = |0| and C^ the number of p-uplets with n numbers. The number 1 will 
be still by convention the intersection of all the elements of Q. The codification 
of di n 63 is given by [1 3] for |e| = 3 and [1 2 4 7] for |e| = 4. And the 
codification of 6^ U 6*3 is given by [1 2 3 4 5 7] for |e| = 3 and [1 2 3 4 6 7 9 12] 
for |e| = 4. 

In order to reduce the complexity, especially using more hardware language 
than Matlab, we could use binary numbers instead of the integer numbers. 

The Smarandache's codification [11], was introduce for the enumeration of 
distinct parts of a Venn diagram. li \Q\ ^ n, < i > denotes the part of 6i with 
no covering with other 6j, i ^ j. < ij > denotes the part of 6i n 6j with no 
covering with other parts of the Venn diagram. So if rt = 2, 611162 = {< 12 >} 
and if n = 3, 01 n 6*2 = {< 12 >, < 123 >}, see the figure 3 for an illustration 
for n — 3. The authors note a problem for n > 10, but if we introduce space in 
the codification we can conserve integers instead of other symbols and we write 
< 1 2 3 > instead of < 123 >. 

On the contrary of the Smarandache's codification, the proposed codification 
gives only one integer number to each part of the Venn diagram. This codifica- 
tion is more complex for the reader then the Smarandache's codification. Indeed, 
the reader can understand directly the Smarandache's codification thanks to the 
mining of the numbers knowing the n: each disjoint part of the Venn diagram is 
seen as an intersection of the elements of Q. More exactly, this is a part of the 
intersections. For example, 61 n 62 is given with the Smarandache's codification 
by {< 12 >} if n = 2 and by {< 12 >, < 123 >} if n == 3. With the codification 
practical codification the same element has also different codification according 
to the number n. For the previous example 61 n 62 is given by [1] if n = 2, and 
by [1 2] if n = 3. 



The proposed codification is more practical for computing union and inter- 
section operations and the DSm cardinahty, because only one integer represent 
one of the distinct parts of the Venn diagram. With the Smarandache's codi- 
fication computing union and intersection operations and the DSm cardinality 
could be very similar than with the practical codification, but adding a routine 
in order to treat the code of one part of the Venn diagram. 




Figure 3: Smarandache's codification for 8 = {Oi, 02,03}- 

Hence, we propose to use the proposed codification to compute union, in- 
tersection and DSm cardinality, and the Smarandache's codification, easier to 
read, to present the results in order to safe eventually a scan of D^. 

3.2 Adding constraints 

With this codification, adding constraints is very simple and can reduce rapidly 
the number of integers. E.g. assume that in a given application we know 
6*1 n 6*3 = {i.e. 6*1 n 6*3 ^ Df), that means that the integers [1 3] for |e| = 3 
and [12 4 7] for \Q\ = 4 do not exist Q. Hence, the codification of 8 with 
the reduced discernment space, noted 8^, is given respectively for |8| =3 and 
|8|=4by: 

8, = {[2 5], [2 4 6], [4 7]} 

and 

8^ = {[3 6 9 12], [3 5 6 8 10 13], [5 8 11 14], [3 5 9 10 11 15]}. 

Generally we have |8| = |0r|, but it is not necessary if a constraint gives 
9i = 0, with 9i g 8. This can happen in dynamic fusion, if one element of the 
discernment space can disappear. 

Thereby, the introduction of the simple constraint 6*1 n 6*3 = in 8, includes 
all the other constraints that follow from it such as the intersection of all the 
elements of 8 is empty. In [15] all the constraints must be given by the user. 



3.3 Codification of the focal elements 

In I?^ , the codification of the focal elements is given from the reduced discern- 
ment space Qr- The codification of an union of two elements of Q is given by the 
concatenation of the codification of the two elements using Or- The codification 
of an intersection of two elements of O is given by the common numbers of the 
codification of the two elements using 8^- In the same way, the codification of 
an union of two focal elements is given by the concatenation of the codification 
of the two focal elements and the codification of an intersection of two focal 
elements is given by the common numbers of the codification of the two focal 
elements. In fact, for union and intersection operations we only consider one 
element as the set of the numbers given in its codification. 

Hence, with the previous example (we assume ^i n 6*3 = 0, with \Q\ = 3 or 
|e| = 4), if the following elements 6*1 n 612, 6'i U 02 and (6*1 n 6*2) U 613 are some 
focal elements, there are coded for |8| = 3 by: 

010 02 = [2], 

01 U 02 = [2 4 5 6], 

(01 n 02) U 03 = [2 4 7], 

and for |6| = 4 by: 

01 n 02 = [3 6], 

6(1 u 02 = [3 5 68 9 10 12 13], 

(01 n 02) U 03 = [3 5 68 11 14]. 

The DSm cardinality Cm (X) of one focal element X is simply given by 
the number of integers in the codification of X. The DSm cardinality of one 
singleton is given by the equation (9), only if there is none constraint on the 
singleton, and inferior otherwise. 

The previous example with the focal element (0i fl 02) U 03 illustrates well 
the easiness to deal with the brackets in one expression. The codification of the 
focal elements can be made with any brackets. 

3.4 Combination 

In order to manage only the focal elements and their associated basic belief 
assignment, we can use a list structure [9, 15, 22]. The intersection and union 
operations between two focal elements coming from two mass functions are made 
as described before. If the intersections between two focal elements is empty the 
associated codification is [ ]. Hence the conjunctive combination rule algorithm 
can be done by the algorithm 1. The disjunctive combination rule algorithm is 
exactly the same by changing n in U. 

Once again, the interest of the codification is for the intersection and union 
operations. Hence in Matlab, we do not need to redefine these operations as in 

[15]. 

For more complicated combination rules such as PCR6, we have generally 
to conserve the intermediate calculus in order to transfer the partial conflict. 
Algorithms for these rules have been proposed in [22], and Matlab codes are 
given in section 5. 



Algorithm 1: Conjunctive rule 



Data: n experts ex: ex[l] . . .ex[n\, ex[i]. focal, ex[i].bba 
Result: Fusion of ex by conjunctive rule: conj 
extmp ^- ex[l]; 
for e = 2 to n do 
comb ^- 0; 

foreach foci in extmp. focal do 
foreach foc2 in ex[e]. focal do 

tm,p ^- extmp. focal (foci) n ex[e].focal{foc2); 

comb. focal ^- tm,p; 

comb.bba ^ extmp.bba{f ocl) x ex[e].bba{foc2); 



Concatenate same focal in comb; 
extmp ^- comb; 
conj <— extm,p; 



3.5 Decision 

As we write before, we can decide with one of the functions given by the equa- 
tions (6), (7), or (8). These functions arc increasing functions. Hence generally 
in 2 , the decision is taken on the elements in O by the maximum of these 
functions. In this case, with the goal to reduce the complexity, we only have 
to calculate these functions on the singletons. However, first, we can provide 
a decision on any element of 2^^ such as in [2] that can be interesting in some 
applications [24], and second, the singletons are not the more precise or inter- 
esting elements on D^ . The figures 4 and 5 show the DSm cardinality Ca^(A), 
\/X € D^ with respectively |6| = 3 and \Q\ — 4. The specificity of the single- 
tons (given by the DSm cardinality) appears at a central position in the set of 
the specificities of the elements in D®. 

Hence, to calculate these decision functions on all the reduced hyper power 
set could be necessary, but the complexity could not be inferior to the complexity 
of Df and that can be a real problem. The more reasonable approach is to 
consider either only the focal elements or a subset of Df on which we calculate 
decision functions. 

3.5.1 Extended weighted approach 

Generally in 2®, the decisions are only made on the singletons [8, 34], and only 

few approaches propose a decision on 2 . In order to provide decision on any 

elements of Df , we can first extend the principle of the proposed approach in 

[2] on Df. This approach is based on the weighting of the plausibility with a 

Bayesian mass function taking into account the cardinality of the elements of 
26 

In a general case, if there is no constraint, the plausibility is not interesting 
because all elements contain the intersection of all the singletons of O. According 
the constraints the plausibility could be applied. 

Hence, we generalize here the weighted approach to Df for every decision 
function fd (plausibility, credibility, pignistic probability, ...). We note fwd the 
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Figure 4: DSm cardinality Cm^X), VX e D® with |e| = 3. 



weighted decision function given for aU X E Df by: 

Ud[X) = md{X)fd{X), 
where nid is a basic behef assignment given by: 



(10) 



(11) 



s is a parameter in [0, 1] allowing a decision from the intersection of all the 
singletons (s = 1) (instead of the singletons in 2^) until the total indecision 
(s = 0). Ax allows the integration of the lack of knowledge on one of the 
elements X in Df. The constant Kd is the normalization factor giving by the 
condition of the equation (4). Thus we decide the element A: 



A = argmax/^d(-^), 

X<^Df 



(12) 



If we only want to decide on whichever focal element of Df^ we only consider 
X G Tm and we decide: 



A = argmax/^<i(X), 



with Jwd given by the equation (10) and: 

/ 1 



md{X) = KdX 



X 



\Cm{X) 
s and Kd are both parameters defined above 
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(13) 
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Figure 5: DSm cardinality Cm(^), VX e Z?® with |e| = 4. 



3.5.2 Decision according to the specificity 

The cardinaUty Cm{X) can be seen as a specificity measure of X. The figures 4 
and 5 show that for a given specificity there is different Icind of elements such as 
singletons, unions of intersections or intersections of unions. The figure 6 shows 
well the central role of the singletons (the DSm cardinality of the singletons 
for |0|=5 is 16), but also that there is many other elements (619) with exactly 
the same cardinality. Hence, it could be interesting to precise the specificity of 
the elements on which we want to decide. This is the role of s in the Appriou 
approach. Here we propose to directly give the wanted specificity or an interval 
of the wanted specificity in order to build the subset of D^ on which we calculate 
decision functions. Thus we decide the element A: 



A 



a.rgmax fd{X), 
xes 



(15) 



where fd is the chosen decision function (credibility, plausibility, pignistic prob- 
ability, ...) and 



S ~ {X e Df;mins < Ca<(X) < maxs} , 



(16) 



with mins and vnaxs respectively the minimum and maximum of the specificity 
of the wanted elements. If mins 7^ maxs, if have to chose a pondered decision 
function for fd such as f^d given by the equation (10). 

However, in order to find all X G S we must scan Df. To avoid to scan all 
D^ , we have to find the cardinality of S, but we can only calculate an upper 
bound of the cardinality, unfortunately never reached. Let define the number 
of elements of the Venn diagram ny- This number is given by: 



ny ^Cm [[J(^i 



(17) 



^4=1 
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Figure 6: Number of elements of D^ for |0| = 5, with the same DSm cardinahty. 

where n is the cardinahty of Qr and 9i e 9^- Recah that the DSm cardinahty 
is simply given by the number of integers of the codification. The upper bound 
of the cardinality of S is given by: 



\S\< E C'^w 



(18) 



where C^ is the number of combinations of s elements among ny ■ Note that 
it also works if mins = for the empty set. 

3.6 Generation of Df 

The generation of Df could have the same complexity than the generation of 
D if there is none constraint given by the user. Today, the complexity of the 
generation of D® is the complexity of the proposed code in [11]. Assume for 
example, the simple constrain 9i Ci 02 = 0. First, the figures 7(a) and 7(b) 
show the DSm cardinality for the elements of Df with |9| = 4 and the previous 
given constraint. On the left figure, the elements are ordered by increasing DSm 
cardinality and on the right figure with the same order than the figure 5. We can 
observe that the cardinality of the elements have naturally decreased and the 
number of non empty elements also. This is more interesting if the cardinality 
of 9 is higher. Figure 8 presents for a given positive DSm cardinality, the 
number of elements of Df for |9| = 5 and with the same constraint 91(102 = ^■ 
Compared to the figure 6, the total number of non empty elements (the integral 
of the curve) is considerably lower. 

Thus, we have to generate Df and not D^ . The generation of D^ (see 
[11] for more details) is based on the generation of monotone boolean functions. 
A monotone boolean function /,„(, is a mapping of {xi, ...,xi,) e {0,1}'' to a 
single binary output such as Vx, x' e {0,1}'', with x =<; x' then fmb{^) < 
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elements of D elements of D® 

(a) Elements are ordered by increasing DSm car- (b) Elements are ordered with the same order 
dinality. than the figure 5. 

Figure 7: DSm cardinality Cm(X), VX G Df with |e| = 4 and 6*1 n 6*2 = 0. 



fmb{^')- Hence, a monotone boolean function is defined by the values of the 
2^ elements {xi, ...,Xb), and there is \D^\ different monotone boolean functions. 
All the values of all these monotone boolean function can be represented by a 
\D''\ X 2^ matrix. If we multiply this matrix by the vector of all the possible 
intersections of the singletons in Q with |0| = 6 (there is 2^ intersections) given 
an union of intersections, we obtain all the elements of D^ . We can also use 
the basis of all the unions of 6 (and obtain the intersections of unions), but 
with our codification the unions are coded with more integer numbers. So, the 
intersection basis is preferable. 

Moreover, if we have some constraints (such as 6*1 fl 02 = 0), some elements 
of the intersection basis can be empty. So we only need to generate a |D*| x nt 
matrix where nt is the number of non empty intersections of elements in Qr- 
For example, with the constraint given in example for |0| = 3, the basis is given 
by: 0, 6*1, 6*2, 6*3, 61 n 63, 6*2 n 63, and there is no 6*1 n 62 and 6*1 n 6*2 n 63. 



Hence, the generation of D^ can run very fast if the basis is small, 
there is some constraints. The Matlab code is given in section 5. 



if 



3.7 Decoding 

Once the decision on one element A of Df is taken, we have to transmit this 
decision to the human operator. Hence we must to decode the element A (given 
by the integer numbers of the codification) in terms of unions and intersections 
of elements of 9. If we know that ^ is in a subset of elements of Df given by 
the operator, we only have to scan this subset. Now, if the decision A comes 
from the focal elements (a priori unknown) or from all the elements of D^ we 
must scan all D^ with possibly high complexity. What we propose here is to 
consider the elements of Df ordering with first the elements most encountered 
in applications. Hence, we first scan the elements of 2® and in the same time 
the intersection basis that we must build for the generation of D^ . Then, only 
if the element is not found we generate D^ and stop the generation when found 
(see the section 5 for more details). 

The Smarandache's codification is an alternative to the decoding because 
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Figure 8: Number of elements of D^ for |0| = 5 and diOO-z ^ ^, with the same 
positive DSm cardinahty. 

user can directly understand it. Hence we can represent the focal element as an 
union of the distinct part of the Venn diagram. The Smarandache's codification 
allows a clear understanding of the different parts of the Venn diagram on the 
contrary than the proposed codification. This representation of the results (for 
the combination or the decision) does not need the generation of D^ . However, 
if we need to generate Df according to the strategy of decision, the decoding 
will give a better display without more generation of Df . 

4 Concluding remarks 

This chapter presents a general belief function framework based on a practical 
codification of the focal elements. First the codification of the elements of the 
Venn diagram gives a codification of O. Then, the eventual constraints are 
integrated giving a reduced discernment space 0,-. From the space 0,-, we 
obtain the codification of the focal elements. Hence, we manipulate elements 
of a reduced hyper power set D^ and not the complete hyper power set D®, 
reducing the complexity according to the kind of given constraints. 

With the practical codification, the step of combination is easily made using 
union and intersection functions. 

The step of decision was particularly studied, because of the diflfculties to 
decide on D^ or Df . An extension of the approach given in [2] in order to give 
the possibility to decide on the unions in 2® was proposed. Another approach 
based on the specificity was proposed in order to simply choose the elements on 
which decide according to their specificity. 

The principal goal of this chapter is to provide practical codes of a general 
belief function framework for the researchers and users needing the belief func- 
tion theory. However, for sake of clarity, all the Matlab codes are not in the 
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listing, but can be provided on demand to the author. The proposed codes 
are not optimized either for Matlab, or in general and can still have bugs. All 
suggestions in order to improve them are welcome. 

5 Matlab codes 

We give and explain here some Matlab codes of the general belief function 
framework^. Note that the proposed codes are not optimized cither for Matlab, 
or in general. 

First the human operator have to describe the problem (see function 1) 
giving the cardinality of O, the list of the focal elements and the corresponding 
bba for each experts, the eventual constraints (' ' if there is no constraint), 
the list of elements on which he want to obtain a decision and the parameters 
corresponding to the choice of combination rule, the choice of decision criterium 
the mode of fusion (static or dynamic) and the display. When the description 
of the problem is made, he just has to call the fuse function 2. 

Function 1 - Command configuration 

7o description of the problem 

CardTheta=4; 'L cardinality of Theta 

°L list of experts with focal elements and associated bba 

expert (l).focal={'l' 'lu3' '3' 'Iu2u3'}; 

expert (l).bba= [0.5421 0.2953 0.0924 0.0702]; 

expert (2). focal={'l' '2' 'lu3' 'Iu2u3'}; 
expert (2). bba= [0.2022 0.6891 0.0084 0.1003]; 

expert (3). focal={'l' '3n4' 'Iu2u3'}; 
expert (3). bba= [0.2022 0.6891 0.1087]; 

constraint={ ' ln2' 'ln3' '2n3'}; °L set of empty elements 
elemDec={ 'F'} ; °L set of decision elements 

I 

7o parameters 

criteriumComb=l ; 7„ combination citerium 
criteriumDec=0; 'L decision criterium 
mode=' static ' ; °L mode of fusion 
display=3; 7, kind of display 

I 

7o fusion 

fuse(expert, constraint ,CardTheta,criteriumComb,criteriumDec, . . . 
mode , elemDec , display) 



^Copyright © 2008 Arnaud Martin. May be used free of charge. Selling without prior 
written consent prohibited. Obtain permission before redistributing. 
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The first step of the fuse function 2 is the coding. The cardinahty of O 
gives the codification of the singletons of G, thanks to the function 3, then 
we add the constraints to O with the function 4 and obtain 0^. With 8^, the 
function 6 calhng the function 5 codes the focal elements of the experts given by 
the human operator. The combination is made by the function 7 in static mode. 
For dynamic fusion, we just consider one expert with the previous combination. 
In this case the order of the experts given by the user can have an important 
signification. The decision step is made with the function 11. The last step 
concern the display and the hard problem of the decoding. Thus, 4 choices are 
possible: no display, the results of the combination only, the results of decision 
only and both results. These displays could take long time according to the 
parameters given by the human operator. Hence, the results of the combination 
could have the complexity of the generation of Df and must be avoid if the 
user does not need it. The complexity of the decision results could also be high 
if the user does not give the exact set of elements on witch decide, or only the 
singletons with 'S' or on 2® with '2T'. In other cases, with luck, the execution 
time can be short thanks to the function 18. 

Function 2 - Fuse function 

function fuse (expert , constraint ,n, cr iter iumComb,cr iter iumDec, mode, elemDec, display) 

y. To fuse experts' opinions 

7. 

7o fuse (expert , constraint , n , criteriumComb , cr iter iumDec , mode , elemDec , display) 

% 

7o Inputs: 

7o expertC = containt the structure of the list of coded focal elements and 

7o corresponding bba for all the experts 

7o constraint = the empty elements 

7o elemDec = list of elements on which we can decide 

7o n = size of the discernment space 

7o criteriumComb = is the combination criterium 

7o criteriumComb=l Smets criterium 

°/t criteriumComb=2 Dempster-Shaf er criterium (normalized) 

7. criteriumComb=3 Yager criterium 

7o criteriumComb=4 disjunctive combination criterium 

7o criteriumComb=5 Florea criterium 

7o criteriumComb=6 PCR6 

7o criteriumComb=7 Mean of the bbas 

7o criteriumComb=8 Dubois criterium (normalized and 

7. disjunctive combination) 

7. criteriumComb=9 Dubois and Prade criterium (mixt combination) 

7. criteriumComb=10 Mixt Combination (Martin and Dsswald criterium) 

7. criteriumComb=ll DPCR (Martin and Osswald criterium) 

7o criteriumComb=12 MDPCR (Martin and Dsswald criterium) 

7o criteriumComb=13 Zhang's rule 

7. 

7. 

7. criteriumDec = is the combination criterium 
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% criteriumDec=0 maximum of the bba 

% criteriumDec=l maximum of the pignistic probability 

% criteriumDec=2 maximum of the credibility 

% criteriumDec=3 maximum of the credibility with reject 

% criteriumDec=4 maximum of the plausibility 

% criteriumDec=5 Appriou criterium 

% criteriumDec=6 DSmP criterium 

% 

7o mode = 'static' or 'dynamic' 

7o elemDec = list of elements on which we can decide, 

y« or A for all, S for singletons only, F for focal elements only, 

% SF for singleton plus focal elements. Cm for given specificity, 

7, 2T for only 2~Theta (DST case) 

°L display = kind of display 

7„ display = for no display, 

7. display = 1 for combination display, 

7o display = 2 for decision display, 

7o display = 3 for both displays, 

7o display = 4 for both displays with Smarandache codification 

7. 

7o Output : 

7. res = containt the structure of the list of focal elements and 

7o corresponding bbas for the combinated experts 

7. 

7. Copyright (c) 2008 Arnaud Martin 

7 Coding 

[Theta,Scod] =codingTheta(n) ; 

ThetaRed=addConstraint (constraint ,Theta) ; 

expertCod=codingExpert (expert , ThetaRed) ; 

7. 

switch nargin 
case 1:5 

mode=' static' ; 
elemDec=ThetaRed; 
display=4; 
case 6 

elemDec=ThetaRed; 
display=4; 
case 7 

elemDec=string2code (elemDec) ; 
display=4; 
end 

7. 

if (display==l) I I (display==2) I I (display==3) 

[DThetar ,D_n] =generationDThetar (ThetaRed) ; 
else 

switch elemDec{l} 
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case {'A'} 

[DThetar ,D_rL] =generationDThetar (ThetaRed) ; 
otherwise 

DThetar. s={[]}; 
DThetar. c={[]}; 
end 
end 

y. 

"/o Combination 

if strcmp(mode , 'static') 

[expertComb] =combinat ion (expertCod, ThetaRed, criteriumComb) ; 
else "/o dynamic case 

nbexp=size (expertCod, 2) ; 

expertTmp(l)=expertCod(l) ; 

for exp=2:nbexp 

expert Tmp (2) =expertCod(exp) ; 

expert Tmp(l)=combination(expertTmp, ThetaRed, criteriumComb) ; 

end 

expertComb=expertTmp(l) ; 
end 

7o Decision 

[decFocElem] =decision(expertComb , ThetaRed, DThetar . c , criteriumDec , elemDec) ; 

7o Display 
switch display 
case 

'no display' 
case 1 

7o Result of the combination 
sFocal=size (expertComb .focal , 2) ; 

focalRec=decodingExpert (expertComb, ThetaRed, DThetar) ; 
f ocal=code2string(f ocalRec) 
for i=l:sFocal 

disp ( [ focal{i}, '=' ,num2str (expertComb. bba(i) ) ] ) 
end 
case 2 

7, Result of the decision 
if isstruct (decFocElem) 

f ocalDec=decodingFocal (decFocElem .focal , elemDec , ThetaRed) ; 
disp( ['decision: ' , code2string(f ocalDec)] ) 
else 

if decFocElem==0 

disp( ['decision: rejected']) 
else 

if decFocElem==-l 

disp( ['decision: cannot be taken']) 
end 
end 
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end 
case 3 

7, Result of the combination 
sFocal=size(expertComb.f ocal,2) ; 

expertDec=decodingExpert (expertComb , ThetaRed , DThetar) ; 
f ocal=code2string(expertDec .focal) 
for i=l:sFocal 

disp ( [ f ocal{i}, '=' ,num2str (expertDec .bba(i) ) ] ) 
end 

7, Result of the decision 
if isstruct (decFocElem) 

focalDec=decodingFocal(decFocElem. focal, elemDec,ThetaJled, DThetar) ; 
disp( ['decision: ' , code2string(f ocalDec)] ) 
else 

if decFocElem==0 

disp( ['decision: rejected']) 
else 

if decFocElem==-l 

disp( ['decision: cannot be taken']) 
end 
end 
end 
case 4 

7o Results with Smarandache codification display 
7o Result of the combination 
sFocal=size (expertComb. focal, 2) ; 
expertDec=cod2ScodExpert (expertComb , Scod) ; 
for i=l:sFocal 

disp ( [expertDec .f ocal{i}, ' = ' ,num2str (expertDec .bba(i) )] ) 
end 

7o Result of the decision 
if isstruct (decFocElem) 

f ocalDec=cod2ScodFocal (decFocElem .focal , Scod) ; 
disp( ['decision: ' ,f ocalDec] ) 
else 

if decFocElem==0 

disp( ['decision: rejected']) 
else 

if decFocElem==-l 

disp( ['decision: cannot be taken']) 
end 
end 
end 

otherwise 

'Accident in fuse: choice of display is uncorrect' 
end 
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5.1 Codification 

The codification is based on the function 3. The order of the integer numbers 
could be different, here the choice is made to number the intersection of all the 
elements with 1 and the smallest integer among the \&\ = n bigger integers 
for the first singleton. In the same time this function give the correspondence 
between the integer numbers of the practical codification and the Smarandache's 
codification. This function 3 is based on the Matlab function nchoosek(tah,k) 
given the array of all the combination of k elements of the vector tah. If the 
length of tah is n, this function return an array of C^ rows and k columns. 

Function 3 - codingTheta function 
function [Theta.Scod] =codingTheta(n) 

7o Code Theta for DSmT framework 

7. 

7o [Theta,Scod] =codingTheta(n) 

7. 

1 Input: 

7o n = cardinality of Theta 

7. 

7o Outputs : 

7o Theta = the liste of coded elements in Theta 

7o Scod = the bijection function between the integer of 

7o the coded elements in Theta and the Smarandache codification 

I 

7 Copyright (c) 2008 Arnaud Martin 

i=2-n-l; 
tablnd= [] ; 
for j=n:-l:l 

tablnd= [tabind j] ; 

Theta{j}=[i] ; 

Scod{i}=[j] ; 
i=i-l; 
end 

i=i+l; 

for card=2:n 

tabPerm=nchoosek (tabind, card) ; 
for j=l:n 

[1, c] =f ind(tabPerm==j) ; 
tabi=i . *ones(l , sized, 1) ) ; 
Theta{j}=[sort(tabi-l') Theta{j}] ; 
for nb=l : sized, 1) 

Scod{i-l(nb)}=[Scod{i-l(nb)} j] ; 
end 
end 
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i=i-size(tabPerm, 1) ; 
end 



The addition of the constraints is made in two steps: first the codification 
of the elements in the hst constraint is made with the function 5, then the 
integer numbers in the codification of the constraints are suppressed from the 
codification of O. The function string2code is just the translation of the brackets 
and union and intersection operators in negative numbers (-3 for '(', -4 for ')', -1 
for 'U' and -2 for 'n') in order to manipulate faster integers than strings. This 
simple function is not provided here. 



Function 4 - addConstraint function 



function [ThetaR] =add.Constraint (constraint ,Theta) 

7o Code ThetaR the reduced form of Theta 

7o taking into account the constraints given by the user 

7. 

7o [ThetaR] =addConstraint (constraint , Theta) 

% 

°L Inputs : 

7o constraint = the list of element considered as constraint 

7. or '2T' to work on 2"Theta 

7o Theta = the description of Theta after coding 

I 

7 Output: 

7 ThetaR = the description of coded Theta after reduction 

7 taking into account the constraints 

7 

7 Copyright (c) 2008 Arnaud Martin 

if strcmp(constraint{l}, '2T') 
n=size (Theta, 2) ; 
nbCons=l ; 
for i=l:n 

for j=i+l:n 

constraint (nbCons)={ [i -2 j]}; 
nbCons=nbCons+l ; 
end 
end 
else 

constraint=string2code (constraint) ; 
end 

constraintC=codingFocal (constraint , Theta) ; 

sConstraint=size(constraintC,2) ; 

unionCons= [] ; 

for i=l : sConstraint 

unionCons=union(unionCons, constraintC{i}) ; 
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end 

sTheta=size(Theta,2) ; 
for i=l:sTheta 

ThetaR{i}=setdif f (Theta{i},unioiiCoiis) ; 
end 



The function 5 simply transforms the hst of focal elements given by the 
user with the codification of 8 to obtain the list of constraints and with 8^ for 
the focal elements of each expert. The function 6 prepares the coding of focal 
elements and return the list of the experts with the coded focal elements. 



Function 5 - codingFocal function 



function [f ocalC] =codingFocal(f ocal,Theta) 

7o Code the focal element for DSmT framework 

7. 

7o [f ocalC] =codingFocal(f ocal,Theta) 

7. 

7o Inputs : 

7o focal = the list of focal element for one expert 

7o Theta = the description of Theta after coding 

I 

7 Output: 

7 focalC = the list of coded focal element for one expert 

7 

7 Copyright (c) 2008 Arnaud Martin 

nbf oc=size (focal , 2) ; 
if nbfoc 

for foc=l:nbfoc 

elemC=treat (focal{foc}, Theta) ; 
f ocalC{f oc}=elemC; 
end 
else 

focalC={[]}; 
end 
end 

77 

function [elemE] =eval(oper ,a,b) 

if oper==-2 

elemE=intersect (a,b) ; 
else 

elemE=union(a,b) ; 
end 
end 
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u 

function [elemC, cmp] =treat (focal, Theta) 

nbelem=size(f ocal,2) ; 

PelemC=0; 

oper=0; 

e=l; 

if nbelem 

while e <= nbelem 
elem=f ocal(e) ; 
switch elem 
case -1 

oper=-l ; 
case -2 

oper=-2; 
case -3 

[elemC,nbe] =treat (focal (e+1 : end) , Theta) : 
e=e+nbe ; 

if oper~=0 & ~isequal(PelemC,0) 

elemC=eval (oper , PelemC , elemC) ; 
oper=0; 

end 

PelemC=eleiiiC; 

case -4 

cmp=e ; 
e=nbelem; 
otherwise 

elemC=Theta{elem} ; 

if oper~=0 & ~isequal (PelemC, 0) 

elemC=eval (oper , PelemC , elemC) ; 
oper=0; 
end 

PelemC=elemC; 
end 
e=e+l; 
end 
else 

elemC= [] ; 
end 
end 



Function 6 - codingExpcrt function 
function [expertC] =codingExpert (expert , Theta) 

7o Code the focal element for DSmT framiework 
I 
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7» [expertC] =codingExpert (expert ,Theta) 

7. 

°L Inputs : 

7» expert = structure containing the list of focal elements for 

7o each expert and the bba corresponding 

7o Theta = the description of Theta after coding 

7. 

7. Output: 

7. expertC = structure containing the list of coded focal element 

7o for each expert and the bba corresponding 

% 

7. Copyright (c) 2008 Arnaud Martin 

nbExp=size (expert ,2) ; 
for exp= 1 : nbExp 

focal=string2code (expert (exp) .focal) ; 

expertC (exp) .focal=codingFocal (focal, Theta) ; 

expertC (exp) .bba=expert (exp) .bba; 
end 
end 



5.2 Combination 

The function 7 proposes many combination rules. Most of them are based on the 
function 8, but for some combination rules we need to keep more information, 
so we use the function 9 for the conjunctive combination. E.g. in the func- 
tion 10 note the simplicity of the code for the PCR6 combination rule. Other 
combination rules' codes are not given here for the sake of clarity. 

Function 7 - combination function 

function [res] =combinat ion ( expertC, ThetaR,cr it erium) 

7o Give the combination of many experts 

7. 

7o [res] =combination(expert .constraint ,n,criterium) 

I 

°L Inputs : 

7o expertC = containt the structure of the list of focal elements 

7o and corresponding bba for all the experts 

7. ThetaR = the coded and reduced discernment space 

7. criterium = is the combination criterium 

7o criterium=l Smets criterium (conjunctive rule in open world) 

7o criterium=2 Dempster-Shaf er criterium (normalized) 

7o (conjunctive rule in closed world) 

7o criterium=3 Yager criterium 

7o criterium=4 disjunctive combination criterium 

7o criterium=5 Florea criterium 

7 criterium=6 PCR6 

7o criterium=7 Mean of the bbas 
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7o criterium=8 Dubois criterium 

7o (normalized and disjunctive combination) 

7o criterium=9 Dubois and Prade criterium (mixt combination) 

°L criterium=10 Mixt Combination (Martin and Osswald criterium) 

°L criterium=ll DPCR (Martin and Osswald criterium) 

7o criterium=12 MDPCR (Martin and Osswald criterium) 

°L criterium=13 Zhang's rule 

7. 

°L Output : 

7o res = containt the structure of the list of focal elements and 

7» corresponding bbas for the combinated experts 

% 

7. Copyright (c) 2008 Arnaud Martin 

switch criterium 
case 1 

7.Smets criterium 
res=conjunctive(expertC) ; 
case 2 

7oDempster-Shaf er criterium (normalized) 
expConj=conjunctive(expertC) ; 
ind=f indeqcell(expConj .focal, [] ) ; 
if "isempty (ind) 

k=expConj .bba(ind) ; 
expConj .bba=expConj .bba/(l-k) ; 
expConj .bba(ind)=0; 
end 

res=expConj ; 
case 3 

7oYager criterium 
expConj=conjunctive(expertC) ; 
ind=findeqcell (expConj .focal, [] ) ; 
if ~isempty (ind) 

k=expConj .bba(ind) ; 
eTheta=ThetaR{l}; 
for i=2:n 

eTheta= [union(eTheta,ThetaR{i})] ; 
end 

indTheta=findeqcell (expConj .f ocal,eTheta) ; 
if ~isempty(indTheta) 

expConj .bba(indTheta)=expConj .bba(indTheta)+k; 
expConj .bba(ind)=0; 
else 

sFocal=size (expConj .focal, 2) ; 
expConj .f ocal(sFocal+l)={eTheta}; 
expConj .bba(sFocal+l)=k; 
expConj .bba(ind)=0; 
end 
end 
res=expConj ; 
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case 4 

Zdisjounctive criterium 

[res] =disjunctive(expertC) ; 
case 5 

7o Florea criterium 

expCoiij=conjunctive(expertC) ; 

expDis=disjunctive(expertC) ; 

ind=f indeqcell(expConj .focal, [] ) ; 
if ~isempty(ind) 

k=expConj .bba(ind) ; 

alp]ia=k/(l-k+k*k); 

beta=(l-k)/(l-k+k*k) ; 

expFlo=expConj ; 
expFlo .bba=beta. *expFlo .bba; 
expFlo . bba ( ind) =0 ; 
nbFocCoiij=size(expConj . focal, 2) ; 
iibFocDis=size (expDis . focal , 2) ; 

expFlo . focal (nbFocConj+1 :nbFocConj+nbFocDis)=expDis . focal; 
expFlo .bba(nbFocConj+l :nbFocConj+nbFocDis)=alpha. *expDis .bba; 

expFlo=reduceExpert (expFlo) ; 

else 

expFlo=expCoiij ; 

end 

res=expFlo; 
case 6 

7. PCR6 

[res] =PCR6 (expertC) ; 
case 7 

7o Means of the bba 

[res] =meanbba(expertC) ; 
case 8 

7o Dubois criterium (normalized and disjunctive combination) 

expDis=disjunctive(expertC) ; 

ind=f indeqcell(expDis . focal, [] ) ; 
if ~isempty(ind) 

k=expDis .bba(ind) ; 
expDis .bba=expDis .bba/(l-k) ; 
expDis .bba(ind)=0; 
end 

res=expDis; 
case 9 

7o Dubois and Prade criterium (mixt combination) 
[res]=DP(expertC) ; 
case 10 

7. Martin and Dssawald criterium (mixt combination) 
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[res] =Mix (expertC) ; 
case 11 

7. DPCR (Martin and Osswald criterimn) 

[res]=DPCR(expertC) ; 
case 12 

7, MDPCR (Martin and Osswald criterium) 

[res] =MDPCR(expertC) ; 
case 13 

7. Zhang's rule 

[res] =Z]iang( expert) 

otherwise 

'Accident: in combination choose of criterium: uncorrect' 
end 



Function 8 - conjunctive function 
function [res] =conjunctive (expert) 

°L Conjunctive Rule 

7. 

7o [res] =conjunctive (expert) 

7. 

7o Inputs : 

7o expert = containt the structures of the list of focal element and 

7o corresponding bba for all the experts 

7. 

7. Output: 

7. res = is the resulting expert (structure of the list of focal 

7o element and corresponding bba) 

I 

7. Copyright (c) 2008 Arnaud Martin 

nbexpert=size (expert , 2) ; 
for i=l:nbexpert 

nbf ocal(i)=size(expert (i) .focal, 2) ; 
nbbba(i)=size(expert (i) .bba, 2) ; 
if nbf ocal(i) ~=nbbba(i) 

'Accident: in conj : the numbers of bba and focal element... 
are different' 
end 
end 

interm=expert (1) ; 
for exp=2:nbexpert 

nbf ocallnterm=size(interm. focal, 2) ; 

i=l; 

comb. f ocal={} ; 

comb . bba= [] ; 
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for f ocl=l :nbf ocallnterm 

for f oc2=l :nbf ocal(exp) 

tmp=intersect (interin.focal{f ocl}, expert (exp) .f ocal{f oc2}) ; 
if isempty (tmp) 

tmp= [] ; 
end 

comb . focal (i) ={tmp} ; 

comb. bba(i)=iiiterm.bba(f ocl) *expert (exp) .bba(f oc2) ; 
i=i+l; 
end 
end 

interm=reduceExpert (comb) ; 
end 
res=interm; 



Function 9 - globalConjunctivc function 
function [res,tablnd] =globalConjunctive (expert) 

7o Conjunctive Rule conserving all the focal elements 

7o during the combination 

7. 

7o [res ,tablnd] =globalConjunctive (expert) 

7. 

1 Input: 

7o expert = containt the structures of the list of focal element and 

7o corresponding bba for all the experts 

I 

°L outputs : 

7o res = is the resulting expert (structure of the list of focal 

7o element and corresponding bba) 

7o tabind = table of the indices given the combination 

I 

7. Copyright (c) 2008 Arnaud Martin 

nbexpert=size (expert , 2) ; 
for i=l:nbexpert 

nbfocal(i)=size (expert (i) . focal, 2) ; 

nbbba(i)=size(expert (i) .bba, 2) ; 

if nbf ocal(i) ~=nbbba(i) 

'Accident: in conj : the numbers of bba and focal element... 
are different' 

end 
end 

interm=expert (1) ; 
tabIndPrev=[l:l:nbfocal(l)] ; 
for exp=2:nbexpert 

nbf ocallnterm=size (interm . focal , 2) ; 
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i=l; 

comb. f ocal={} ; 

comb . bba= [] ; 

tablnd=[] ; 

for f ocl=l :nbf ocallnterm 

for f oc2=l :nbf ocal(exp) 

tinp=iiitersect (interm.focal{f ocl}, expert (exp) .f ocal{f oc2}) ; 
tabliid=[tabliid [tablndPrevC : ,focl) ;foc2]] ; 
if isempty(tmp) 

tmp= [] ; 
end 

comb . focal (i) ={tmp} ; 

comb.bba(i)=interm.bba(f ocl)*expert (exp) .bba(foc2) ; 
i=i+l; 
end 
end 

tabIndPrev=tabInd ; 
interm=comb; 
end 
res=interm; 



Function 10 - PCR6 function 

function [res] =PCR6 (expert) 

7o PCR6 combination rule 

7. 

7, [res] =PCR6 (expert) 

% 

1 Input: 

7o expert = containt the structures of the list of focal element and 

7o corresponding bba for all the experts 

I 

7 Output: 

7 res = is the resulting expert (structure of the list of focal 

7. element and corresponding bba) 

7. 

7 Reference: A. Martin and C. Osswald, ''A new generalization of the 

7 proportional conflict redistribution rule stable in terms of decision, ' ' 

7 Applications and Advances of DSmT for Information Fusion, Book 2, 

7 American Research Press Rehoboth, F. Smarandache Eind J. Dezert, 

7 pp. 69-88 2006. 

7 

7 Copyright (c) 2008 Arnaud Martin 

[expertConj ,tablnd] =globalConjunctive (expert) ; 

ind=f indeqcell(expertConj .focal, [] ) ; 
nbexp=size(tablnd, 1) ; 
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if ~isempty(ind) 

expertConj .bba(ind)=0; 
slnd=size(ind,2) ; 
for i=l:slnd 
P=l; 
S=0; 
for exp=l:nbexp 

bbaexp=expert (exp) .bba(tabliid(exp, ind(i)) ) ; 
P=P*bbaexp; 
S=S+bbaexp; 
end 
for exp=l:iibexp 

expertConj .f ocal(end+l)=expert (exp) .f ocal(tabInd(exp, ind(i) )) ; 
expertConj .bba(end+l)=expert (exp) .bba(tablnd(exp, ind(i)) )*P/S; 
end 
end 
end 
res=reduceExpert (expertConj) ; 

5.3 Decision 

The function 11 gives the decision on the expert focal element list for the cor- 
responding bba with one of the chosen criterium and on the elements given by 
the user for the decision. Note that the choices 'A' and 'Cm' for the variable 
elemDec could take a long time because it need the generation of Df. This 
function can call one of the decision functions 13, 14, 15, 16. If any decision is 
possible on the chosen elements given by elemDec, the function return -1. In 
case of reject element, te function return 0. 

Function 11 - decision function 

function [decFocElem] =decision(expert ,Theta, criterium, elemDec) 

7o Give the decision for one expert 

7. 

7o [decFocElem] =decision (expert ,Theta, criterium) 

7. 

7o Inputs : 

7o expert = containt the structure of the list of focal elements and 

7o corresponding bba for all the experts 

7o Theta = list of coded (and reduced with constraint) of the 

7o elements of the discernement space 

7o criterium = is the combination criterium 

7o criterium=0 maximum of the bba 

7o criterium=l maximum of the pignistic probability 

7o criterium=2 maximum of the credibility 

7o criterium=3 maximum of the credibility with reject 

7o criterium=4 maximum of the plausibility 

7o criterium=5 DSmP criterium 

7o criterium=6 Appriou criterium 
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7o criterium=7 Credibility on DTheta criterium 

y, criterium=8 pignistic on DTheta criterium 

7, elemDec = list of elements on which we can decide, 

7o or A for all, S for singletons only, F for focal elements only, 

7o SF for singleton plus focal elements. Cm for given specificity, 

7. 2T for only 2~Theta (DST case) 

7. 

7 Output: 

7 decFocElem = the retained focal element, in case of rejet, -1 

7o if the decision cannot be taken on elemDec 

I 

7. Copyright (c) 2008 Arnaud Martin 



type=l; 

switch elemDec{l} 
case 'S' 

type=0; 

elemDecC=Theta; 
expertDec=expert ; 
case 'F' 

elemDecC=expert . focal ; 
expertDec=expert ; 
case 'SF' 

expertDec=expert ; 
n=size(Theta,2) ; 
for i=l:n 

expertDec . f ocal{end+l}=Theta{i} ; 
expertDec .bba(end+l)=0; 
end 

expertDec=reduceExpert (expertDec) ; 
elemDecC=expertDec . focal ; 
case 'Cm' 

sElem=size (elemDec, 2) ; 
switch sElem 
case 2 

minSpe=str2num(elemDec{2}) ; 
maxSpe=minSpe ; 
case 3 

minSpe=str2num(elemDec{2}) ; 
maxSpe=str2num(elemDec{3}) ; 
otherwise 

'Accident in decision: with the option Cm for . 
elemDec give the specifity of decision element 
(eventually the minimum and the maximum of the 
desired specificity' 
pause 
end 



elemDecC=f indFocal(Theta,minSpe,maxSpe) ; 
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expertDec . f ocal=elemDecC ; 

expertDec .bba=zeros(l , size(elemDecC,2) ) ; 
for f oc=l : size(expert . f ocal,2) 

ind=findeqcell(elemDecC, expert .f ocal{f oc}) ; 
if "isempty (ind) 

expertDec .bba(ind)=expert .bba(f oc) ; 
else 

expertDec . bba ( ind) =0 ; 
end 
end 

case '2T' 
type=0; 

natoms=size(Theta,2) ; 
expertDec .f ocal(l)={ [] }; 
indFoc=findeqcell (expert .f ocal,{ [] }) ; 
if isempty (indFoc) 

expertDec .bba(l)=0; 
else 

expertDec .bba(l)=expert .bba(indFoc) ; 
end 

step =2; 
for i=l:natoms 

expertDec . focal (step) =codingFocal({ [i] },Theta) ; 

indFoc=f indeqcell (expert . focal , expertDec . f ocal{step}) ; 
if isempty (indFoc) 

expertDec .bba(step)=0; 
else 

expertDec .bba(step)=expert .bba(indFoc) ; 
end 

step=step+l ; 

indatom=step; 

for step2=2 : indatom-2 

expertDec . focal (step) ={ [union(expertDec .f ocal{step2}, . 

expertDec . f ocal{indatom-l})] }; 

indFoc=f indeqcell (expert .focal, expertDec .f ocal{step}) ; 
if isempty (indFoc) 

expertDec .bba(step)=0; 
else 

expertDec .bba(step)=expert .bba(indFoc) ; 
end 

step=step+l ; 
end 
end 
elemDecC=expertDec . focal ; 
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case 'A' 

elemDecC=generationDThetar (Theta) ; 
elemDecC=reduceFocal(elemDecC) ; 
expertDec . f ocal=elemDecC ; 
expertDec .bba=zeros(l , size(elemDecC,2) ) ; 
for f oc=l : sizeCexpert . f ocal,2) 

expertDec .bba(findeqcell(elemDecC, expert .f ocal{f oc}) ) 
=expert .bba(f oc) ; 
end 
otherwise 
type=0; 

eleiiiDec=string2code(elemDec) ; 
elemDecC=codiiigFocal(elemDec, Theta) ; 
expertDec=expert ; 
iibEleiiiDec=size(elemDecC,2) ; 
for f oc=l :nbElemDec 

if ~isElem(elemDecC{f oc}, expertDec . focal) 
expertDec .f ocal{eiid+l}=elemDecC{f oc}; 
expertDec .bba(end+l)=0; 
end 
end 
end 

I 

nbFocal=size (expertDec . focal , 2) ; 
switch criterium 

case 

7o maximum of the bba 
nbFocal=size (expertDec .focal, 2) ; 
nbElem=0 ; 
for f oc=l :nbFocal 

ind=findeqcell(elemDecC, expertDec .f ocal{f oc}) ; 

if ~isempty(ind) 

bba(ind)=expertDec .bba(f oc) ; 

end 
end 

[bbaMax , indMax] =max (bba) ; 
if bbaMax~=0 

decFocElem . bba=bbaMax ; 

decFocElem . f ocal={elemDecC{indMax}} ; 
else 

decFocElem=-l ; 
end 
case 1 

7o maximum of the pignistic probability 
[BetP] =pignistic (expertDec) ; 
decFocElem=MaxFoc (BetP , elemDecC , type) ; 
case 2 

7, maximum of the credibility 
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[Bel] =credibility(expertDec) ; 
decFocElem=MaxFoc (Bel , elemDecC , type) ; 
case 3 

7, maximum of the credibility with reject 
[Bel] =credibility(expertDec) ; 

TabSing=[] ; 

f ocTheta= [] ; 

for i=l : size(Theta,2) 

f ocTheta=uiiioii(f ocTheta,Theta{i}) ; 
end 

for f oc=l :nbFocal 

if isElem(Bel . f ocal{f oc}, elemDecC) 

Tabs ing= [Tabs ing [foe ; Bel.BeKf oc)] ] ; 

end 
end 

[BelMax,indMax]=max(TabSing(2, :)) ; 
if BelMax~=0 

f ocMax=Bel . f ocal{TabSing(l , indMax) } ; 

f ocComplementary=setdif f (f ocTheta,f ocMax) ; 
if isempty (f ocComplementary) 

f ocComplementary= [] ; 
end 

ind=f indeqcell(Bel .focal,f ocComplementary) ; 

if BelMax < Bel.Bel(ind) 

7, if ind is empty this is always false 

decFocElem=0; °h That means that we reject 
else 

if isempty (ind) 

decFocElem=0; % That means that we reject 
else 
decFocElem.f ocal={Bel .f ocal{TabSing(l , indMax)}}; 
decFocElem . Bel=BelMax ; 
end 
end 
else 

decFocElem=-l ; 7. That means that we reject 
end 
case 4 

7. maximum of the plausibility 
[PI] =plausibility(expertDec) ; 
decFocElem=MaxFoc (PI , elemDecC , type) ; 
case 5 

7o DSmP criterium 

epsilon=0 .00001 ; °h can allows problem 

[DSmP] =DSmPep (expertDec , epsilon) ; 
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decFocElem=MaxFoc (DSmP , elemDecC , type) ; 
case 6 

7, Appriou criterium 

[PI] =plausibility(expertDec) ; 

lajiibda=l ; 

r=0.5; 

bm=BayesianMass (expertDec , lambda , r ) ; 

Newbba=Pl . PI . *bm . bba; 

7o normalization 

Newbba=Newbba/sum(Newbba) ; 

f uncDec . f ocal=Pl . focal ; 

funcDec .bba=Newbba; 

decFocElem=MaxFoc (f uncDec , elemDecC , type) ; 
case 7 

7o Credibility on DTheta criterium 

[Bel] =credibility (expertDec) ; 

lainbda=l ; 

r=0.5; 

bm=BayesianMass (expertDec , lambda , r ) ; 

Newbba=Bel . Bel . *bm . bba ; 

7o normalization 

Newbba=Newbba/sum(Newbba) ; 

f uncDec . f ocal=Bel .focal ; 

f uncDec .bba=Newbba; 

decFocElem=MaxFoc (f uncDec , elemDecC , type) ; 
case 8 

7o pignistic on DTheta criterium 

[BetP] =pignistic (expertDec) ; 

lainbda=l ; 

r=0.5; 

bm=BayesianMass (expertDec , lambda , r ) ; 

Newbba=BetP . BetP . *bm . bba ; 

7. normalization 

Newbba=Newbba/sum(Newbba) ; 

f uncDec . f ocal=BetP . focal ; 

f uncDec .bba=Newbba; 

decFocElem=MaxFoc (f uncDec , elemDecC , type) ; 
otherwise 

'Accident: in decision choose of criterium: uncorrect' 
end 
end 



VL 

function [bool] =isElem(f ocal, listFocal) 

7o The g oal of this function is to return a boolean on the test focal in 

7o listFocal 

7. 

7o [bool] =isElem(f ocal, listFocal) 
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7. 

7o Inputs : 

7o focal = one focal element (matrix) 

°L listFocal = the list of elements in Theta (all different) 

y. 

°L Output : 

7o bool = boolean: true if focal is in listFocal, elsewhere false 

y. 

7. Copyright (c) 2008 Arnaud Martin 

n=size (listFocal , 2) ; 
bool=f alse ; 
for i=l:n 

if isequal(listFocal{i}, focal) 
bool=true; 
break; 
end 
end 
end 

7.7. 

function [decFocElem] =MaxFoc (f uncDec , elemDecC , type) 

f ieldN=f ieldncmies(funcDec) ; 

switch fieldN{2} 
case 'BetP' 

f uncDec . bba=f uncDec . BetP ; 
case 'Bel' 

f uncDec . bba=f uncDec . Bel ; 
case 'PI' 

f uncDec . bba=f uncDec . PI ; 
case 'DSmP' 

f uncDec . bba=f uncDec . DSmP ; 

end 

if type 

[f uncMax, indMax]=max(f uncDec .bba) ; 
FocMax={f uncDec . f ocal{indMax}} ; 
else 

nbFocal=size (f uncDec . focal , 2) ; 

TabSing= [] ; 

for f oc=l :nbFocal 

if isElem(f uncDec .focal{foc}, elemDecC) 

Tabs ing= [Tabs ing [foe ; funcDec .bba(f oc)] ] ; 
end 
end 

[funcMax,indMax]=max(TabSing(2, :)) ; 
FocMax={f uncDec . f ocal{TabSing(l , indMax)}} ; 
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end 

if funcMax~=0 

decFocElem . focal =FocMax ; 
switch fieldN{2} 
case 'BetP' 

decFocElem . BetP=f uncMax ; 
case 'Bel' 

decFocElem . Bel=f uncMax ; 
case 'PI' 

decFocElem. Pl=f uncMax; 
case 'DSmP' 

decFocElem . DSmP=f uncMax ; 
end 
else 

decFocElem=-l ; 
end 

end 



Function 12 - findFocal function 



function [elemDecC] =f indFocal(Theta,minSpe,maxSpe) 

7o Find the element of DTheta with the minium of specifity minSpe 

7o and the maximum maxSpe 

7. 

7o [elemDecC] =f indFocal(Theta,minSpe,maxSpe) 

7. 

7 Input: 

7 Theta = list of coded (and eventually reduced with constraint) of 

7 the elements of the discernment space 

7 minSpe = minimum of the wanted specificity 

7 minSpe = maximum of the wanted specificity 

7 

7 Output: 

7 elemDec = list of elements on which we want to decide with the 

7 minimum of specifity minSpe and the maximum maxSpe 

7 

7 Copyright (c) 2008 Arnaud Martin 

elemDecC{l}=[] ; 

n=size (Theta, 2) ; 
ThetaSet= [] ; 
for i=l:n 

ThetaSet=union(ThetaSet ,Theta{i}) ; 
end 
for s=minSpe : maxSpe 
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tabs=nchoosek(ThetaSet ,s) ; 

elemDecC(end+l : end+size(tabs , 1) )=num2cell(tabs,2) ' ; 
end 

elemDecC=elemDecC(2: end) ; 



Function 13 - pignistic function 
function [BetP] =pignistic (expert) 

7o Generalized Pignistic Transformation 

7. 

7o [BetP] =pignistic (expert) 

7. 

7. Input: 

7« expert = containt the structures of the list of focal element and 

7 corresponding bba for all the experts 

7 expert. focal = list of focal elements 

7 expert. bba = matrix of bba 

7 

7 Output : 

7 BetP = containt the structure of the list of focal element and 

7 the matrix of the plausibility corresponding 

7 BetP. focal = list of focal elements 

7 BetP. BetP = matrix of the pignistic transformation 

7 Comment : 1- the code of the focal elements must inculde 
% the constraints 

% 2- The pignistic is given only on the elements 
% in the list of focal of expert (the 

% bba can be 0) 

7 

7 Copyright (c) 2008 Arnaud Martin 

nbFocal=size (expert . focal , 2) ; 

BetP . f ocal=expert . focal ; 
BetP.BetP=zeros(l,nbFocal) ; 

for f ocA=l :nbFocal 

for f ocB=l :nbFocal 

f ocI=intersect (expert . f ocal{f ocA} , expert . f ocal{f ocB}) ; 
if "isempty (f ocl) 

BetP.BetP(focA)=BetP.BetP(focA)+size(focI,2)/. . . 
size (expert . f ocal{f ocB},2)*expert .bba(f ocB) ; 
else 

if isequal (expert . f ocal{f ocB} , [] ) 
% for the empty set : 
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°/o cardinality (empty set) /cardinality (empty set)=l, 
7o so we add the bba 

BetP . BetP (f ocA) =BetP . Bet? (f ocA) +expert . bba(f ocB) ; 
end 
end 
end 
end 



Function 14 - credibility function 
function [Bel] =credibility (expert) 

7o Credibility function 

7. 

7o [Bel] =credibility (expert) 

7. 

1 Input: 

7« expert = containt the structures of the list of focal element and 

7« corresponding bba for all the experts 

7« expert. focal = list of focal elements 

7« expert. bba = matrix of bba 

7. 

7. Output: 

7o Bel = containt the structure of the list of focal element and 

7o the matrix of the credibility corresponding 

7o Bel. focal = list of focal elements 

7o Bel. Bel = matrix of the credibility 

7o Comment : 1- the code of the focal elements must inculde 

7o the constraints 

7. 2- The credibility is given only on the elements 

7. in the list of focal of expert (the 

7. bba can be 0) 

7. 

7. Copyright (c) 2008 Arnaud Martin 

nbFocal=size (expert . focal , 2) ; 

Bel . f ocal=expert . focal ; 
Bel .Bel=zeros(l ,nbFocal) ; 

for f ocA=l :nbFocal 

for f ocB=l :nbFocal 

indMem=ismember (expert .focal{focB}, expert .f ocal{f ocA}) ; 

if sum (indMem)==size (expert .f ocal{f ocB},2) 

Bel.Bel(focA)=Bel.Bel(focA)+expert.bba(focB) ; 
else 

if isequal (expert . f ocal{f ocB} , [] ) 
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7o the empty set is include to all the focal elements 

Bel.Bel(focA)=Bel.Bel(focA)+expert.bba(focB) ; 
end 
end 
end 
end 



Function 15 - plausibility function 
function [PI] =plausibility (expert) 

7o Plausibility function 

7. 

7o [PI] =plausibility (expert) 

% 

1 Input: 

7o expert = containt the structures of the list of focal element and 

7o corresponding bba for all the experts 

7o expert. focal = list of focal elements 

7o expert. bba = matrix of bba 

7 

7 Output: 

7 PI = containt the structure of the list of focal element ajid 

7 the matrix of the plausibility corresponding 

7 PI. focal = list of focal elements 

7 PI. PI = matrix of the plausibility 

7 Comment : 1- the code of the focal elements must inculde 

7 the constraints 

7o 2- The plausibility is given only on the elements 

7, in the list of focal of expert (the 

% bba can be 0) 

7 

7 Copyright (c) 2008 Arnaud Martin 



nbFocal=size (expert . focal , 2) ; 
PI . f ocal=expert .focal ; 
Pl.Pl=zeros(l,nbFocal) ; 

for f ocA=l :nbFocal 

for f ocB=l :nbFocal 

focI=intersect (expert . f ocal{f ocA} , expert . f ocal{f ocB}) ; 
if ~isempty(f ocl) 

PI . PI (f ocA) =P1 . PI (f ocA) +expert . bba (f ocB) ; 
else 

if isequal (expert . f ocal{f ocB} ,[])... 

&& isequal (expert . f ocal{f ocA} , [] ) 
% for the empty set we keep the bba for the PI 
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PI . PI (f ocA) =P1 . PI (f ocA) +expert . bba (f ocB) ; 
end 
end 
end 
end 



Function 16 - DSmPcp function 
function [DSmP] =DSmPep (expert , epsilon) 

y. DSmP Transformation 

7. 

7, [DSmP] =DSmPep (expert , epsilon) 

% 

°L Inputs : 

7o expert = containt the structures of the list of focal element and 

7o corresponding bba for all the experts 

7o expert. focal = list of focal elements 

7o expert. bba = matrix of bba 

7o epsilon = epsilon coefficient 

7 

7 Output: 

7 DSmPep = containt the structure of the list of focal element and 

7 the matrix of the plausibility corresponding 

7 DSmPep. focal = list of focal elements 

7 DSmPep. DSmP = matrix of the pignistic transformation 

7 

7 Reference: Dezert & Smarandache , 

7 ''A new probbilistic transformation of belief mass assignment'', 

7 fusion 2008, Cologne, Germany. 

7 

7 Copyright (c) 2008 Arnaud Martin 



nbFocal=size (expert . focal , 2) ; 

DSmP . f ocal=expert . focal ; 
DSmP.DSmP=zeros(l,nbFocal) ; 

for f ocA=l :nbFocal 

for f ocB=l :nbFocal 

f ocI=intersect (expert .focal{focA}, expert .f ocal{f ocB}) ; 

sumbbaFocB=0 ; 

sFocB=size (expert .f ocal{f ocB},2) ; 

for elB=l:sFocB 

ind=findeqcell (expert .focal, expert .f ocal{f ocB}(elB)) ; 
if "isempty (ind) 

sumbbaFocB=sumbbaFocB+expert .bba(ind) ; 
end 
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end 

if "isempty (f ocl) 
siimbbaFocI=0; 
sFocI=size(f ocl ,2) ; 
for elB=l:sFocI 

ind=findeqcell (expert .f ocal,f ocI(elB)) ; 
if ~isempty(ind) 

sumbbaFocI=sumbbcLFocI+expert .bba(iiid) ; 
end 
end 

DSmP.DSmP(focA)=DSmP.DSmP(focA)+expert.bba(focB) . . . 
* (sumbbaFocI+epsilon*sFocI)/(sumbbaFocB+epsilon*sFocB) ; 
end 
end 
end 



5.4 Decoding and generation of Df 

For the displays, wc must decode the focal elements and/or the final decision. 
The function 17 decodes the focal elements in the structure expert that contain 
normally only one expert. This function calls the function 18 that really does 
the decoding for the user. This function is based on the generation of Df given 
by the function 21 that a is modified and adapted code from [11]. To generate 
-D® we first must create the intersection basis. Hence in the function 18 we use 
a loop of 2^^ in order to generate the basis and in the same time to scan the 
power set 2® and also the elements of the intersection basis. These two basis 
(intersection and union) are in fact concatenated during the construction, so 
we scan also some elements such intersections of previous unions and unions 
of previous intersections. This generated set of elements does not cover Df. 
When all the searching focal elements (that can be only one decision element) 
are found, we stop the function and avoid to generate all Df. Hence if the 
searching elements are not all found after this loop, we begin to generate Df 
and stop when all elements arc found. So, with luck, that can be fast. 

We can avoid to generate Df for only the display if we use the Sniarandachc's 
codification. The function 19 transforms the used code of the focal elements in 
the structure expert in the Smarandache's code, easer to understand by reading. 
This function calls the function 20 that really docs the transformation. The focal 
elements are directly in string for the display. 



Function 17 - decodingExpert function 



function [expertDecod] =decodingExpert (expert ,Theta,DTheta) 

y. The goal of this function is to decode the focal elements in expert 

7. 

7» [expertDecod] =decodingExpert (expert ,Theta) 

t 

°L Inputs: 

°L expert = containt the structure of the list of focal elements after 

7o combination and corresponding bba for all the experts (generally use 
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°L for only one after combination) 

7o Theta = list of coded (and reduced with constraint) of the elements of 

y. the discernement space 

7o DTheta = list of coded (and reduced with constraint) of the elements of 

% DTheta 

% 

7, Output: 

°L expertDecod = containt the structure of the list of decoded (for human) 

7o focal elements and corresponding bba for all the experts 

7. 

7 Copyright (c) 2008 Arnaud Martin 

nbExp=size (expert ,2) ; 
for exp= 1 : nbExp 

f ocal=expert (exp) . focal ; 

expertDecod (exp) . f ocal=decodingFocal(f ocal,{' A'},Theta,DTheta) ; 

expertDecod (exp) .bba=expert (exp) .bba; 
end 

end 



Function 18 - decodingFocal function 
function [f ocalDecod] =decodingFocal(focal,elemDec, Theta, DTheta) 

7 The goal of this function is to decode the focal elements 

7 

7 [f ocalDecod] =decodingFocal (focal , elemDec , Theta) 

7 

7 Inputs : 

7 expert = containt the structure of the list of focal elements after 

7 combination and corresponding bba for all the experts 

7 elemDec = the description of the subset of uncoded elements 

7 for decision 

7 Theta = list of coded (and reduced with constraint) of the 

7 elements of the discernement space 

7 DTheta = list of coded (and reduced with constraint) of the 

7 elements of DTheta, eventually empty if not necessary 

7 Output: 

7 focalDecod = containt the list of decoded (for human) focal elements 

7 

7 Copyright (c) 2008 Arnaud Martin 

switch elemDec{l} 

case {'F' , 'A', 'SF', 'Cm'} 

opt=l; 
case 'S' 

opt=0; 

elemDecC=Theta ; 
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for i=l : size(Theta,2) 

elemDec(i)={ [i] }; 
end 
case '2T' 
opt=0; 

natoms=size(Theta,2) ; 
elemDecC (!)={[]}; 
elemDec(l)={[]}; 
step =2; 
for i=l:iiatoms 

elemDecC (step) =codingFocal({ [i] },Theta) ; 

elemDec (step) ={ [i] } ; 

step=step+l ; 

indatom=step; 

for step2=2 : iiidatoiii-2 

elemDec (step) ={ [elemDec{step2} -1 elemDec{indatom-l}] }; 
elemDecC (step) ={ [union (elemDecC{step2} , elemDecC{indatom-l}) ] } ; 
step=step+l ; 

end 
end 

otherwise 
opt=0 ; 

elemDecN=string2code (elemDec) ; 
elemDecC=codingFocal(elemDecN,Theta) ; 
end 

if "opt 

sFoc=size(f ocal,2) ; 
for foc=l:sFoc 

[ind] =findeqcell (elemDecC, focal{foc}) ; 
if isempty(ind) 

'Accident in decodingFocal : elemDec does not be 2T' 
pause 
else 

f ocalDecod(f oc)=elemDec(ind) ; 
end 
end 
else 

focalDecod=cell (size (focal) ) ; 
cmp=0 ; 

sFocal=size(f ocal,2) ; 
sDTheta=size(DTheta. c,2) ; 
i=l; 

while KsDTheta && cmp<sFocal 
DThetai=DTheta . c{i} ; 
indeq=findeqcell (focal, DThetai) ; 
if "isempty (indeq) 
cmp=cmp+l ; 
f ocalDecod(indeq)=DTheta. s(i) ; 
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end 
i=i+l; 
end 
end 



Function 19 - cod2ScodExpert function 



function [expertDecod] =cod2ScodExpert (expert , Scod) 

7o The goal of this function is to code the focal elements in 

7o expert with the Smarandache ' s codification from the practical 

7o codification in order to display the expert 

7. 

7o [expertDecod] =cod2ScodExpert (expert , Scod) 

% 

°L Inputs : 

7o expert = containt the structure of the list of focal elements after 

7o combination and corresponding bba for all the experts (generally use 

7o for only one after combination) 

7o Scod = list of distinct part of the Venn diagram coded with the 

7o Smarandache ' s codification 

7 Output: 

7 expertDecod = containt the structure of the list of decoded (for human) 

7 focal elements and corresponding bba for all the experts 

7 

7 Copyright (c) 2008 Arnaud Martin 

nbExp=size (expert ,2) ; 
for exp= 1 : nbExp 

f ocal=expert (exp) . focal ; 

expertDecod(exp) . focal=cod2ScodFocal (focal, Scod) ; 

expertDecod (exp) .bba=expert (exp) .bba; 
end 

end 



Function 20 - cod2ScodFocal function 



function [f ocalDecod] =cod2ScodFocal (focal, Scod) 

7 The goal of this function is to code the focal elements with the 

7 Smarandache ' s codification from the practical codification in order to 

7 display the focal elements 

7 

7 [f ocalDecod] =cod2ScodFocal (focal, Scod) 

7 

7 Inputs: 
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7» expert = containt the structure of the list of focal elements after 

7o combination and corresponding bba for all the experts 

% Scod = list of distinct part of the Venn diagram coded with the 

7o SmarEindache ' s codification 

7. Output: 

7o focalDecod = containt the list of decoded (for human) focal elements 

% 

7 Copyright (c) 2008 Arnaud Martin 

sFocal=size(f ocal,2) ; 
for foc=l:sFocal 

sElem=size(f ocal{f oc},2) ; 
if sElem==0 

f ocalDecod{f oc}= ' {} ' ; 
else 

ch='{'; 

ch=strcat (ch, '<' ) ; 

ch=strcat (ch,num2str(Scod{f ocal{f oc}(l)}) ) ; 

ch=strcat (ch, '>' ) ; 

for elem=2:sElem 

ch=strcat (ch, ',<'); 

ch=strcat (ch,num2str (Scod{f ocal{f oc}(elem)}) ) ; 
ch=strcat (ch, '>' ) ; 
end 

f ocalDecod{f oc}=strcat (ch, '}' ) ; 
end 
end 



Function 21 - generationDThetar function 
function [DTheta] =generationDThetar (Theta) 

7 Generation of DThetar: modified and adapted code from 

7 Dezert & Smarandache Chapter 2 DSmT book 7 Vol 1 to generate DTeta 

7 

7 [DTheta] =generationDThetar (Theta) 

7 

7 Input: 

7 Theta = list of coded (and eventually reduced with constraint) of 

7 the elements of the discernment space 

% 

7 Output: 

7 DTheta = list of coded (and eventually reduced with constraint in 

7 this case some elements can be the same) of the elements 

% of the DTheta 

7 

7 Copyright (c) 2008 Arnaud Martin 

n=size (Theta, 2) ; 
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step =1; 
for i=l:n 

basetmp(step)={ [Theta{i}] }; 

step=step+l ; 

indatom=step; 

for step2=l : indatom-2 

basetmp(step)={intersect (basetmp{indatom-l},basetmp{step2})}; 
step=step+l ; 

end 
end 

sBaseTmp=size(basetmp,2) ; 
step=l ; 
for i=l : sBaseTmp 

if "isempty (basetmp{i}) 

base (step) =basetmp(i) ; 
step=step+l ; 

end 
end 

sBase=size (base , 2) ; 
DTheta{l}=[] ; 
step=l ; 
nbC=2; 
stop=0; 

D_nl =[0 ; 1]; 
sDnl=2; 
for nn=l:n 

D_n =[ ] ; 

cf irst=l+(nn==n) ; 
for i =l:sDnl 

Li=D_nl(i, :) ; 
sLi=size(Li ,2) ; 

if (2*sLi>sBase)&& (Li(sLi-(sBase-sLi) )==!) 
stop=l ; 
break 
end 

for j=i:sDnl 

Lj=D_nl(j,:); 

if (and(Li,Lj)==Li)&(or(Li,Lj)==Lj) 
D_n=[D_n ; Li Lj ] ; 

if size(D_n, l)>step 
step=step+l ; 
DTheta{step}=[] ; 
for c=cfirst:nbC 
if D_n(end,c) 

if isempty (DTheta{step}) 

DTheta{step}=base{sBase+c-nbC} ; 
else 
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DTheta{step}=union(DTheta{step},base{sBase+c-nbC}) ; 
end 
end 
end 
end 
end 
end 
end 
if stop 

break 
end 

D_nl=D_n; 

sDnl=size (D_nl , 1) ; 
nbC=2*size (D_nl , 2) ; 
end 
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